Electronic structure of the molecule based magnet Cu PM(NOs)2 (H 2 0)2 



o 
o 



in 

(N 



CO 



K. Doll 

Institut fur Mathematische Physik, TU Braunschweig, 38106 Braunschweig, Germany 
Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569 Stuttgart, Germany 

A. U. B. Wolter 

Institut fur Physik der Kondensierten Materie, TU Braunschweig, 38106 Braunschweig, Germany 
Hahn-Meitner-Institut at BESSY, Albert- Einstein- Strafie 15, 12489 Berlin, Germany 

H.-H. Klauss 

Institut fiir Physik der Kondensierten Materie, TU Braunschweig, 38106 Braunschweig, Germany 

(Dated: February 1, 2008) 

We present density functional calculations on the molecule based S = | antiferromagnetic chain 
compound Cu PM(NOs)2 (H 2 0) 2 ; PM = pyrimidine. The properties of the ferro- and antiferro- 
magnetic state are investigated at the level of the local density approximation and with the hybrid 
functional B3LYP. Spin density maps illustrate the exchange path via the pyrimidine molecule 
which mediates the magnetism in the one-dimensional chain. The computed exchange coupling is 
antiferromagnetic and in reasonable agreement with the experiment. It is suggested that the an- 
tiferromagnetic coupling is due to the possibility of stronger derealization of the charges on the 
nitrogen atoms, compared to the ferromagnetic case. In addition, computed isotropic and anisotropic 
hyperfine interaction parameters are compared with recent NMR experiments. 

PACS numbers: 75.30.Et, 75.50.Ee, 75.50.Xx 
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I. INTRODUCTION 

Molecule based magnetic materials have been the sub- 
ject of intense research, with the target of designing 
new magnetic materials For example, novel quantum 
phenomena such as quantum tunneling of magnetization 
open up possible future applications in quantum com- 
puting and data storage @,H, 0, H[ . Other interesting re- 
search directions are spin state transitions or tuning 
of the magnetic coupling and the design of ferromagnets 
(see, e.g. 0, 9)- This has led to an interdisciplinary 
effort in physics and chemistry, and by experimentalists 
and theoreticians. 

Supramolecular complexes of transition metals with 
organic ligands can also be used to synthesize low- 
dimensional magnets. The organic ligands constitute 
magnetic superexchange pathways with a strength of the 
order of 1 to 100 Kelvin. Since pyrimidine and simi- 
lar heterocycles (pyrazine, pyridine) are often found as 
magnetic exchan ge m ediating molecules in metal-organic 
magnets (e.g. [9ITT0I. [TTt]^. it is important to study in de- 
tail the electronic structure and magnetic exchange mech- 
anism. 

Cu PM(NOs)2 (H 2 0)2 is a molecule based magnet 
which can be considered as a one dimensional spin chain 
(see figure [1} . This complex was synthesized a few years 
ago [I4I. It has been studied by magnetic suscepti- 
bility, specific heat and electron spin resonance mea- 
surements |13| . More recently, high field magnetiza- 
tion studies 11411 and 13 C NMR measurements were per- 
formed [l5|, [HI. It can be theoretically described as a 
S = h antiferromagnetic Heisenberg chain with an ex- 
change coupling of J=36 K [TU, [lj, QjJ, and with an 




FIG. 1: A cutout of the one dimensional chain of Cu 
PM(N0 3 ) 2 (H 2 0) 2 . 



additional Dzyaloshinskii-Moriya interaction and a stag- 
gered g-tensor 0, [l8[ . This model had also been used for 
other one-dimensional spin chain systems such as copper 
benzoate [HHil or CuCl 2 • 2(dimethylsulfoxide) 

Earlier, a molecular orbital study based on the ex- 
tended Hiickel approach had been performed to gain an 
insight in the origin of the magnetic interaction and to 
study the magnetic pathway [221 ] . The extended Hiickel 
method can be viewed as a first step in the hierarchy of 
ab initio calculations. 

In this article, we present a density functional study 
of this system, in order to obtain results based on first 



2 



principles calculations, without using experimental data 
(apart from the positions of the nuclei). The target is to 
get an understanding of the charge- and spin distribution 
by an analysis of spin density maps and a calculation of 
the individual magnetic moments. By computing the en- 
ergy difference between ferro- and antiferromagnet, the 
exchange coupling J can be extracted. In addition, NMR 
parameters such as the isotropic and anisotropic hyper- 
fine interaction parameters are computed and compared 
with recent experimental values. The aim is to get an un- 
derstanding of the counter-intuitive experimental result, 
that an atom with a relatively large distance to the mag- 
netic ion has a larger isotropic shift than an atom closer 
to this magnetic ion. The density functional approach 
allows to obtain all these properties on equal footing. 

II. METHOD 

The calculations were done with the code CRYS- 
TAL2003 [23| . This code employs a local basis set made 
of Gaussian type functions. For Cu, a [5s4p2ef] basis set 
[24]], for O a [As3p] basis set[25| (with outermost sp ex- 
ponents of 0.5 and 0.191), for N as [3s2pld], for C a 
[3s2pld], and for H a [2slp] basis set was chosen; the 
latter three basis sets were as in (26|. Full potential, 
all electron density functional calculations with the local 
density approximation (LDA) and with the hybrid func- 
tional B3LYP were performed. These calculations were 
done for the ferro- and antiferromagnetic state, where 
the resultant solution of the Kohn-Sham equations is an 
eigenstate of S z , but not of S 2 . The energy difference was 
therefore fitted to an Ising model, in order to estimate 
the exchange coupling J. From the computed spin den- 
sity, the isotropic and the anisotropic hyperfine coupling 
parameters are extracted. The charge and spin of the in- 
dividual atoms are obtained via the Mulliken population 
analysis. 

III. RESULTS 

A. Charge and spin densities 

In tabic [H the Mulliken populations of the ferromag- 
netic solution are displayed. Copper carries a charge 
of ~ +1.6, i.e. less than a formal charge of +2. Con- 
sequently, the total spin is ~ 0.7, which indicates that 
the spin is delocalized to the neighboring atoms. Con- 
cerning the pyrimidine ring, we notice that the nitrogen 
atoms are negatively (-0.7) and the carbon atoms posi- 
tively charged, so that the ring as a whole is positively 
charged (0.3). The charge on NO3 is ~ -0.9, and H 2 is 
approximatively neutral. The largest spin on the pyrim- 
idine ring (~ 0.1) is located on the nitrogen atoms of the 
pyrimidine ring which are neighbors to the copper ions. 
Comparing LDA and B3LYP, we note that the LDA so- 
lution gives a slightly more delocalized picture. This is 



TABLE I: Results from the Mulliken population analysis for 
the ferromagnetic solution. 

B3LYP LDA 



atom 


charge 


spin 


charge 


spin 


Cu 


1.6 


0.7 


1.5 


0.6 


N 


-0.7 


0.09 


-0.6 


0.128 
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H bonded to CI 


0.03 


0.003 


0.04 


0.003 


C2 


0.09 


0.01 


0.10 


0.02 


H bonded to C2 


0.02 


0.002 


0.03 


0.003 


C3 


0.4 


-0.009 


0.4 


-0.005 


H bonded to C3 


0.000 


0.002 


-0.001 


0.003 


=>• Pyrimidine 


0.3 


0.2 


0.4 


0.3 


NO3 


-0.9 





-0.8 


0.01 


H2O 


-0.1 


0.05 


-0.1 


0.06 



consistent with previous findings, e.g. [13, [H, [2{| |3(| 
where it was shown that LDA overemphasizes dereal- 
ization. The spin in the pyrimidine ring is alternating up 
and down, consistent with the idea of a spin polarization 
mechanism. Looking at the individual sites, for CI and 
C3 essentially the p-orbitals carry the spin, whereas for 
C2 the C s orbital carries a little more spin. 

In table HH the corresponding results for the spin of the 
antiferromagnetic solution are displayed. The charges 
are virtually identical to the charges of the ferromagnetic 
solution and thus not displayed. The total spin is similar 
to the ferromagnetic case for the copper atom, and for 
the nitrogen atoms of the pyrimidine ring (apart from 
the sign, obviously). The spin is zero for the CI and C2 
sites due to the symmetry. For the C3 site, there is in the 
case of the LDA a very small spin, parallel to the spin of 
the nearest copper, in contrast to the ferromagnetic case, 
where the spin is antiparallcl. 

These results are visualized in the spin density plots 
for ferromagnetic (figure [2]) and antiferromagnetic (figure 
[3]) spin density, at the B3LYP level. It is obvious that 
neighboring copper and nitrogen atoms always have par- 
allel spin. In the ferromagnetic case, the adjacent carbon 
atoms (C3) have a spin which is antiparallel to the nitro- 
gen atoms: antiparallcl spin between neighboring atoms 
reduces the Pauli repulsion. This allows a more diffuse 
charge distribution and thus reduces the energy. Finally, 
the carbon atom C2 has again an antiparallel spin with 
respect to the neighboring carbon atoms (C3), in order 
to reduce the Pauli repulsion. In the antiferromagnetic 
case, the arguments hold similarly, but additionally, the 
cancellation of negative and positive spin density must be 
taken into account which results in a total spin of zero 
or nearly zero at the carbon sites. 

Comparing the ferro- and antiferromagnetic solutions, 
the antiferromagnetic spin density allows the charge den- 
sities of the two nitrogen atoms to stronger interpenetrate 
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TABLE II: Results from the Mulliken population analysis for 
the antiferromagnetic solution. 





B3LYP 


LDA 


atom 


spin 


spin 


Cu a,b 


± 0.7 


± 0.6 


N a,b 


± 0.08 


± 0.10 


CI 


0.000 


0.000 


H bonded to CI 


0.000 


0.000 


C2 


0.000 


0.000 


H bonded to C2 


0.000 


0.000 


C3 a,b 


0.000 


± 0.003 


H bonded to C3a,b 


± 0.001 


± 0.002 


=> Pyrimidine 








NO3 bonded to Cu a,b 





± 0.01 


H2O attached to Cu a,b 


± 0.05 


± 0.06 



TABLE III: Fermi contact coupling, at the level of LDA and 
B3LYP, versus the experimental value, in 1 / a§ , with the Bohr 
radius ao. 

site LDA B3LYP experiment [16] 
CI 0.008 0.003 0.0045 
C2 0.008 0.007 0.034 
C3 0.001 -0.0006 -0.006 




in the center of the ring, which is suppressed in the fer- 
romagnetic case because of the Pauli principle. This is 
illustrated in figure [5] where it is shown that the charge 
density in the ring is slightly higher for the antiferromag- 
netic solution. As a whole, this results in a stronger der- 
ealization of the nitrogen charge in the antiferromagnetic 
case and thus reduces the total energy, which supports 
the antiferromagnetic coupling observed experimentally. 



FIG. 2: Spin density of the ferromagnetic solution, in atomic 
units (l/a , with the Bohr radius ao), at the B3LYP level. 
Full contour lines represent positive spin density, dashed lines 
negative spin density, and the dashed-dotted line represent 
zero spin density. The lines have a distance of 0.00025/aQ. 



B. Magnetic hyperfine interaction 

From the spin densities, it is possible to compute the 
isotropic Fermi contact and the dipolar contribution to 
the anisotropic hyperfine coupling. The isotropic part is 
given by the spin density at the carbon nuclei, and can 
be compared with the values obtained from 13 C NMR 

0. As a magnetic field is applied in the experiment, we 
therefore have to use the Fermi couplings obtained with 
the ferromagnetic solution. The data are displayed in 
Table \TU\ It should be mentioned that computing Fermi 
contact couplings accurately is a notorious problem al- 
ready for molecules [3~fl | . This is even more difficult here, 
as the spin density must be evaluated at the position of 
atoms which are far away from the magnetic copper ions, 

1. e. transferred hyperfine fields. The B3LYP approach re- 



produces the signs of the spin densities at the different 
atoms properly: the spin density at the carbon nucleus is 
positive for CI and C2, and negative for C3. Note that 
the total spin is negative for CI, whereas the spin density 
at the nucleus is positive and even in rough quantitative 
agreement with the experimental value. For the C3 site, 
we find a small negative spin density at the nucleus, in 
qualitative agreement with the experiment. Finally, for 
the C2 site, a positive Fermi contact coupling is found 
which is the largest of the values computed. Again, this 
is in qualitative agreement with the experiment. Note 
that the value is largest at this site which has the largest 
distance to the copper atoms; one might rather expect 
to find larger Fermi couplings for atoms with shorter dis- 
tances. From the spin density plots, we find an inter- 
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FIG. 3: Spin density of the antiferromagnetic solution, in 
atomic units (1/aj)), at the B3LYP level. Full contour lines 
represent positive spin density, dashed lines negative spin den- 
sity, and the dashed-dotted line represent zero spin density. 
The lines have a distance of 0.00025/aQ. 



FIG. 4: The difference of the charge density of the ferro- 
magnetic and the antiferromagnetic solution, in atomic units 
(|e|/ao), at the B3LYP level. Full contour lines represent re- 
gions where the density is higher in the ferromagnetic case, 
dashed lines represent regions where the density is higher 
in the antiferromagnetic case. The dashed-dotted line cor- 
responds to zero. The lines have a distance of 0.00001|e|/ao . 



pretation and explanation: first, in the pyrimidine ring, 
neighboring atoms have antiparallel spin, so that we find 
alternating up and down spin. As the nitrogen spin is 
very large, the CI atom has a relatively small spin density 
because the negative CI spin is partially compensated by 
the positive nitrogen spin which is spatially extended to- 
wards the CI site. In addition, the carbon spin resides 
essentially in the p orbital which has a node at the nu- 
cleus and thus does not contribute to the spin density. 
In contrast, at the C2 site, a spin parallel to the nitrogen 
spin is obtained and mediated via both adjacent carbon 
atoms. The carbon s orbital, which has a non-vanishing 
spin density at the carbon nucleus, contributes slightly 
more to the spin for the C2 site. Thus, we find a rel- 
atively large isotropic hyperfine coupling although this 



atom has the largest distance to the copper atoms car- 
rying the majority of the spin. Finally, the carbon atom 
at C3 has relatively small spin with opposite sign; here 
again the negative spin of the carbon is compensated by 
the neighboring positive spin of the nitrogen and C2. 

It becomes also apparent that B3LYP fits better to the 
experimental values than LDA does; the different sign for 
the site C3 can only be confirmed at the B3LYP level. 

In a next step, the components of the anisotropic hy- 
perfine tensor are computed. These are computed as the 
expectation value of the operator 
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where ip u are the Gaussian type basis functions, P^^g 
is the density matrix for the difference of up and down 
spin, g are lattice vectors, and Ta is the distance r-A, 
with A being the position of the nucleus for which the 
anisotropic dipole hyperfine tensor is computed. The re- 
sults of these calculations are presented in table HVl A 
comparison is made with the results from a point lat- 
tice dipole moment where 90 % of the spin was allocated 
at the copper site, and 5 % at each of the two nitro- 
gen atoms of the pyrimidine ring. This was found to 
be the best fit to the measured NMR data [l6j]. Such 
a spin distribution is similar to the one obtained in the 
present work by a first principles density functional calcu- 
lation. Comparing the dipolar tensors, we find a reason- 
able agreement between the three approaches used (point 
dipolar model, LDA, B3LYP). This is also demonstrated 
in figure where the computed anisotropic hyperfine in- 
teraction tensor was used, together with experimentally 
determined NMR chemical shift and isotropic hyperfine 
interaction. It can be seen that the data is reasonably 
well fitted, but also the dependence on the functional 
becomes obvious; and B3LYP fits the data better. 

It should be mentioned that the results for the hy- 
perfine interaction were obtained with the basis set as 
described in section ITT] As this property depends on the 
spin density of the nucleus, additional calculations were 
performed with e.g. a much larger set of tighter carbon s 
and p basis functions to describe the electron density in 
the region of the nucleus better. We found that the re- 
sults were essentially stable with respect to various basis 
sets employed. 



C. Magnetic exchange coupling 

The magnetic exchange coupling is usually obtained by 
computing the energy difference between two magnetic 
states, and fitting it to a model Hamiltonian. This means 
that the energies of at least two magnetic states have 
to be computed, in order to extract J from the energy 
difference. 

In the case of molecules, the energy can be com- 
puted with accurate wave function based methods such 
as the multi-reference configuration interaction scheme, 
multi-reference perturbation theory and related methods 
[H S3, [H, HIH SI]. The advantage is that the cor- 
relation treatment is well controlled, and by taking into 
account more and more determinants, a fairly accurate 
calculation of the exchange coupling is possible. In addi- 
tion, the wave function can be constructed as an eigen- 
state of S 2 . On the other hand, such calculations are 
very demanding in terms of memory requirement, and 
are usually limited to few magnetic centers (two in most 
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FIG. 5: A comparison of experimental data for the CI 
(squares), C2 (circles) and C3 (triangles) sites and a fit us- 
ing experimental values for NMR chemical shift and isotropic 
shift, and computed values (LDA, B3LYP) for the anisotropic 
shift. 



cases). In addition, it may be necessary to truncate large 
molecules and create some embedding. An alternative is 
to use methods such as unrestricted Hartree-Fock theory 
or density functional theory and apply them to molecules. 

In the case of solids, a wave-function based treatment 
of the periodic solid is usually prohibitive. Thus, a possi- 
ble way is to use schemes such as Hartree-Fock theory or 
density functional theory which can be applied to solids. 
This approach can be performed whenever the unit cell 
is not too large. 

The main downside of Hartree-Fock theory is that it 
gives a too localized picture and usually strongly under- 
estimates the exchange coupling, often by a factor of 3: 
e^g. for NiO and MnO [37|, KXF 3 (X=MnTfe,Co,Ni) 
[II HI] or KCuF 3 , K 2 CuF 4 and Sr 2 Cu0 2 Cl 2 ]M,\M- 

The local density approximation, on the other hand, 
strongly overestimates the exchange coupling, because 
the charges are too delocalized and thus exchange inte- 
grals are too large. This overestimation is often as large 
as a factor of 5: e.g. for NiO [27|, KCuF 3 , K 2 CuF 4 
and Sr 2 Cu0 2 Cl 2 [Z(|. Gradient corrections only slightly 
change this and again, an overestimation was observed, 
e.g. [H, [4l|. The situation is more difficult in sys- 
tems where various couplings are important, e.g. [42j . 
The hybrid functional B3LYP was initially designed for 
molecules, but has become very popular in solid state 
physics because the band gaps obtained are in surpris- 
ingly good agreement with experiment (43|. It interpo- 
lates between Hartree-Fock theory and density functional 
theory, and is now also frequently used for the calcula- 
tion of superexchange coupling constants, where it over- 
estimates exchange couplings by a factor of the order of 
< 2, e.g. NiO H, KCuF 3 , K 2 CuF 4 and Sr 2 Cu0 2 Cl 2 
g3, El, or La 2 Cu0 4 , La 2 Ni0 4 , KNiF 3 , NiF 2 , MnF 2 , 
KMnF 3 Q, or FeCl 2 (PM) 2 and NiCl 2 (PM) 2 0. Fi- 
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nally, it should be mentioned that there are exceptions 
to these simple rules of thumb, especially in cases where 
the two magnetic centers and the bridging atom(s) show 
a strong deviation from an 180° angle and approach 90°, 
i.e. strongly deviate from being on a straight line, e.g. 
Ji in MnF 2 and NiF 2 [13. li^. or the ferric wheel-like 
molecule as in '2'.). [3(j. In these cases, the coupling is 
small according to the Goodenough-Kanamori rules [l| 
and a calculation becomes more difficult. 

A broken symmetry approach and subsequent spin pro- 
jection was suggested as a way of obtaining eigenstates of 
S 2 in the case of molecules j46|. Further suggestions to 
deal with magnetic states were the spin-restricted open 
shell Kohn-Sham (ROKS) [47[ and the spin-restricted 
ensemble- referenced Kohn-Sham method (REKS) [4g| . 
There is an ongoing discussion about the validity of the 
various approaches, see e.g. [SiJ HH [EH- Very re- 
cently, Ruiz et al suggested that one problem was that 
the self interaction error was taken into account twice 
when spin projection was applied together with a self in- 
teraction correction [53| . This was however challenged 
[54j , and it was argued that there was no firm theoretical 
basis for this argument. 

In the case of solids, the only computationally feasible 
way is to use a broken symmetry approach (without spin- 
projection). The solution of the Hartree-Fock or Kohn- 
Sham equations is thus in general not an eigenstate of S 2 , 
but only of S z , and the spatial symmetry is broken. The 
data can be fitted to an Ising model, and should rather 
not be fitted to the Heisenberg model. This approach 
was actually recommended as a 'simple yet elegant way 
out of this problem' (5f| (where 'this problem' refers to 
the problem described in the preceding paragraph). 

A different way of treating solids would be to use some 
embedded cluster scheme which again allows to apply 
the same quantum chemical methods as in the molecular 
case. However, the truncation is not obvious and poses 
again an approximation. 

The issue of using configuration interaction or spin- 
projection was also discussed in the context of quantum 
dots: systems with few electrons are considered, and 
a model Hamiltonian is chosen where parameters such 
as the effective mass and the dielectric constant are ex- 
tracted from the experiment. This allows to construct 
the wave function on the level of configuration inter- 
action and as an eigenstate of S 2 , and the importance 
of doing so was discussed for this class of systems, e.g. 
[H,[53,l5l,[5§|. However, these systems are very different 
from the one considered here: in the present work, local- 
ized spins are considered, and the orbital occupancy of 
the magnetic ions is essentially determined by the crystal 
field. The ground state is thus often better described by 
a single reference wave function, compared to the case of 
molecules, where often a multi-reference wave function is 
necessary. In the case of quantum dots, the situation is 
different: besides examples where the local density ap- 
proximation works surprisingly well, e.g. [60], there are 
other situations where a description by a single deter- 



minant may be poor, e.g. in the case of large magnetic 
fields or in double dots [53] , and it becomes necessary to 
use configuration interaction schemes. 

As we are interested in a uniform description of proper- 
ties such as the spin and charge density and NMR param- 
eters, we therefore evaluated these properties and the ex- 
change coupling at the same level of theory. The strength 
of the magnetic exchange interaction J is thus computed 
by fitting the energy difference between ferro- and anti- 
ferromagnet to an Ising model: H = ^ JSiSi+i- The 
data obtained for J are a prerequisite required if one was 
interested in quantum tunneling; the calculation of the 
anisotropy would be a further step (which requires spin- 
orbit coupling). 

There are two copper atoms per cell, and thus two cou- 
plings of the size J. The energy difference between ferro- 
and antiferromagnet is thus AE = Efm — Eaf = 2zS 2 J 
where z is the number of couplings per cell, i.e. 2 in this 
case. For S = \, we obtain thus AE — J. The computed 
couplings are displayed in table [V] At the B3LYP level, a 
value of 76 K is obtained, at the LDA level a value of 603 
K. The B3LYP value is by a factor of 2 too large, com- 
pared to the experimental value of J=36 K fl3l Il4, [ItJ • 
The LDA value is even larger, because the LDA density 
is much more delocalized. Such overestimations of com- 
puted exchange couplings are typical for the functionals 
employed, as was mentioned above. The B3LYP density 
is more localized and thus a value closer to the experi- 
ment is obtained. These findings are consistent with the 
Mulliken charges in tables|T]and|TTl where a more covalent 
picture was obtained with the LDA. 
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TABLE IV: Dipolar tensors at the carbon sites: the first value corresponds to the level of a localized dipole model, the second 
to LDA and the third one to B3LYP, in 1/aj). The components are given using a cartesian coordinate system. 

CI 

0.01; 0.007; 0.006 0;0;0 0.022; 0.036; 0.033 

0;0;0 -0.011; -0.007; -0.005 0; 0; 

0.022; 0.036; 0.033 0; 0; 0.0004; -0.0006; -0.001 

C2 

0.0005; 0.0005; 0.001 0;0;0 0.0007; 0.01; 0.003 

0;0;0 0.002; 0.0002; -0.001 0; 0; 

0.0007; 0.01; 0.003 0;0;0 0.002; -0.0007; 0.00001 

C3 

-0.004; -0.01; -0.008 ±0.004; ±0.007; ±0.005 0.002; 0.01; 0.01 

±0.004; ±0.007; ±0.005 0.01; 0.02; 0.02 ±0.002; ±0.006; ±0.004 

0.002; 0.01; 0.01 ±0.002; ±0.006; ±0.004 -0.006; -0.01; -0.01 



TABLE V: Energy difference between ferromagnetic and an- 
tiferromagnetic solution and exchange coupling J. 

functional AE (E h ) J (E h ) J (cV) J (K) 
B3LYP 0.00024 0.00024 0.0065 76 
LDA 0.0019 0.0019 0.052 603 



IV. CONCLUSION 

Density functional calculations on the molecule based 
magnet Cu PM(NOs)2 (H20)2 were performed, using the 
local density approximation and the hybrid functional 
B3LYP. The exchange path via the pyrimidine ring was 
analyzed with spin density plots and with a Mulliken spin 
population. The calculations prove a spin transfer from 
the Cu atom to the adjacent nitrogen atoms as had been 
deduced from the NMR experiments. The spin of the 
nitrogen atoms in the pyrimidine ring is parallel to the 
copper spin in all cases. In the ferromagnetic case, the 
spin on the pyrimidine ring is alternating. The carbon 
atoms CI and C3 have a spin essentially in the p orbitals, 
and on the C2 site also a slightly larger spin is found in 
the carbon s orbital. In the antiferromagnetic case, the 
spin is virtually zero on all the carbon atoms. 

In the case of the antiferromagnet, the charges of the 
two nitrogen atoms of the pyrimidine ring can stronger 
interpenetrate. This delocalization reduces the energy 
and explains why antifcrromagnctism is observed. In ad- 
dition, the isotropic and anisotropic hyperfine interaction 
parameters were computed. For the isotropic parame- 
ters, a qualitative agreement could be observed with the 
B3LYP functional. Especially, the experimental result 



that the isotropic hyperfine coupling is largest at the C2 
site, which has the largest distance to the magnetic ion, 
could be confirmed. It is suggested that this is due to the 
relatively large contribution to the spin density from the 
s orbital for the C2 site. In the case of the anisotropic 
dipolar hyperfine tensor, a good agreement with the ex- 
perimental data was found for both functionals, where 
again B3LYP performed better. Finally, the exchange 
coupling was computed via the energy difference between 
ferro- and antiferromagnetic state. A reasonable agree- 
ment was found at the B3LYP level, whereas the local 
density approximation results in by far too large values 
of J due to an enhanced delocalization, which is a well 
known problem of exchange couplings computed with the 
LDA. 
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